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Abstract 

Surface strain rate is best observed by fitting a strain-rate ellipsoid 
to the measured movement of a stake network, or other collection of surface 
features, using a least-squares procedure. Error of the resulting fit 
varies as (LAt/rT)”*, where L is the stake separation. At is the time period 
between initial and final stake survey, and n is the number of stakes in the 
network. This relation suggests that if n is sufficiently high, 
the traditional practice of re-visiting stake-network sites on successive 
field seasons may be replaced by a less costly single-year operation. A 
demonstration using Ross Ice Shelf data shows that reasonably accurate 
measurements can be obtained from 12 stakes after only 4 days of deformation. 
The least-squares procedure may also aid airborne photogrametric 
surveys in that reducing the time interval between survey and re-survey 
could permit better surface-feature recognition. 



INTRODUCTION 


Horizontal strain rate, vertical vorticity and horizontal velocity 
observed at the surface of large ice sheets and ice shelves constitute 
field data essential to the study of large-scale ice-flow dynamics and 
mass balance. These data are traditionally acquired by implanting stake 
networks and surveying their deformation over a known time interval [Nye, 
1959; Zumberg and others, 1960; Thomas and others, in press; Drew and 
Whillans, in press]. In practice, a one-year time interval is required 
to allow relative stake displacements in excess of measurement resolution. 
This one-year period necessitates the deployment of survey parties on two 
successive field seasons. This paper presents a technique, based on 
prior work by Nye [1957, 1959], that reduces the time interval between 
initial and final survey and eliminates costly survey party re-deployment. 

Typical stake networks used to measure the general flow of the Ross 
Ice Shelf and its grounded margins, for example, consist of four individual 
stakes arranged in either a rosette or a rectangular pattern [Thomas and 
others, in press]. These networks provide a degree of 

data redundancy because only three stakes are required to measure horizontal 
strain rates and the vertical component of vorticity. In practice, no 
significance is assigned to apparent deformation gradients within the 
network, so this redundancy is generally treated by averaging the results 
of the four triangular sub-arrays [Thomas and others, in press]. 

This study examines the consequence of expanded data redundancy within 
compact stake networks to determine whether the time interval between 
survey and re-survey can be shortened to less than one field season. 

Figure (1) shows an example of such a stake network. The results 
demonstrate that greater redundancy can shorten field operations and 
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suppress aliasing errors that may occur when small-scale strain-rate 
fluctuations are superimposed on the average large-scale deformation. 

The stake-network design features that achieve these advantages may 
additionally be useful for certain remote sensing programs that employ 
large numbers of natural surface features which can become obscured from 
one year to the next. 

Sources of measurement error: 

The horizontal strain-rate tensor components E ',. (i * 1, 2; j s 1,2) and vertical 
vorticity component are defined by the following expressions: 



where u = (uj, u^, u^) is the ice velocity and subscripts 1, 2, and 3 refer 

to a orthogonal coordinate system (x, y, z) having unit vectors e , e and e 

x y z 

such that e is perpendicular to the geoid and e and e define the horizontal 

x y 

plane. 

Measurement of stake displacements by conventional surveying equipment 
establishes the velocity gradients required to solve Equations (1) and (2) 
for e . . and w . Sufficient time must be allowed between initial and final 
survey so that the displacements will exceed measurement error and any 
unnatural stake disturbances caused by wind or sunlight. If, for example, 

3 

aL = 0.02 m is the stake position uncertainty, L = 1.5 x 10 m is the stake 
separation, and e = 1 x 10”^ see"* is the strain rate scale (typical for 
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Ross Ice Shelf conditions), then a time Interval of 

.* _ 100 aL 

(3) at = — - 0.4 years 

L c 

is required between initial and final stake survey to achieve a displacement 

measurement accuracy of 1%. This time span is generally too long to be fit 

within a single field season; thus a second field season, incurring the 

duplicate cost of field party re-deployment, is required. 

An alternative to allowing one year between survey and re-survey is to 

plant more stakes in each network. From elementary statistics, the 

1 

measurement error should reduce roughly as 7= where n is the number of 
stakes (or where n is the number of independent measurements of a single 
stake) [Mendenhall and Scheaffer, 1973]. It is thus conceivable that, 
by choosing n sufficiently high. At can be reduced to several days or weeks. 

A second potential advantage gained by using stake arrays with large 
numbers of stakes is the reduction of aliasing errors caused by small-scale 
strain rate and vorticity fluctuations superimposed on the large-scale deformation. 
Defining e^. a as the strain rate averaged over a suitably chosen large- 
scale distance and e. . as the fluctuation around this average, the total 

* vJ 

strain rate e. . may be defined by 

* J 


(4) 




e 


v 


ij 


Uniaxial strain-rate data from the Ross Ice Shelf, shown in Figure (2), 

demonstrates that e.. v is not necessarily smaller than 1 . . a and can vary 

* J 1 J 

over a length scale commensurate with the actual 1-2 km stake-separations 
generally used on ice shelves [Thomas and others, in press]. 
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Many field programs are designed to measure e^ 8 rather than c^j V ; 
however. It may not be possible to select a stake placement that will 
average out the effects of the unknown c^ v . Aliasing errors will 
consequently arise In the measured value of e^ 8 . The reduction of 
aliasing errors can be achieved through either Increasing the stake 
separation or Increasing the number of stakes to allow. In effect, 
averaging of individual aliasing errors. Planting stakes with large 
separations may be ruled out by equipment limitations, so Increasing the 
number of stakes may again provide the best alternative. 

Strain-rate estimation from an over-detenained data set 

An n-leg generalization of the 3-leg rosette stake network used by 
Thomas and others [In press] on the Ross Ice Shelf will be examined to 
determine the advantages of redundancy and demonstrate the data-analysls 
technique. Each outlying stake of this rosette defines a position vector (x., y.) 
i = l,n that will change by (6x., 6y ^ ) i = l,n as a result of horizontal 
ice-shelf strain in the time At,. , 1 * l,n between initial survey and re-survey. 
For purposes of this analysis, (i) the net strain is infinitesimal and (ii) 
the strain-rate variations over the length scale spanned by the stake 
network are considered noise. Rosette design procedures other than those 
described here should be adopted if either of these conditions are not met. 

The set of n observed stake displacements are converted into the four 

3U . 

horizontal velocity component gradients — - i,j = 1,2 by solving the 

3x j 

following linear equation: 



(5) 


z 
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/ 3u x 8u x 3u v \ T 

Here s Is the column vector 1 gjj-, •gy-, -gj*, J representing 


the true 

local average of the horizontal velocity gradients, d Is the column vector 
^«x, fix, fix n fiy, fiy, fiy 


«x, fix 2 fix n fiy, «y 2 - 

7T“* 7t^» • ••» — * 7*~* •••» ^t~ J that re P resents the measured stake 


atj* at 2 


At n atj» at 2 


n 


displacements resulting from the true ice strain s, z is the column vector 


1° 1 °_Z_ 

\*i' A V 


'1 *2 u n °n+l V? *2 n 

7t7» atl» ItT’ At, 7t~ J that represents the disturbance of the 

1 c nit / 


true stake displacements resulting from Inaccurate measurement and from 


actual strain rate variation over the rosette, and A is the 2n x 4 stake-position 
matrix given by 


( 6 ) 


x i y i 0 0 

• ft* 

x y_ 0 0 
n n 

0 0 Xj y 1 


The noise vector z is assumed to be a random variable characterized by the 
following 2n x 2n covarience matrix <z z T >: 
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<zz T > 


( 7 ) 




A V 


n+1 

Atj2 


2n 

A V 


In the above representation of < z_ z^>, the error in measuring one stake is 
assumed* for simplicity, to be uncorrelated with those of other stakes. The 
uncertainty of one stake's displacement, however, may be greater than another's. 
In more complex situations, errors in measuring different stakes may be 
correlated. In this circumstance, the data pre-weighting matrix W defined 
in the following section should be chosen to diagonalize < z zj>. 

The matrix A contains all the Information regarding the number and 
positioning of the rosette stakes. In general, A will have a greater number 
of rows than columns because it maps a given element s of the "velocity 
gradient vector space" having a dimensionality of 4 Into an element d of 
the "data vector space" having a dimensionality equal to twice the number of 
stakes. Note, however, that not all possible elements of the data space are 
accessible through A from an element of the velocity gradient space. These 
inaccessible data elements constitute the measurement of noise and provide 
no information useful for calculating The objective of the data analysis 
technique presented here is to select an inverse of A that discards all 
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Inaccessible data elements. 

Pre-weighting and scaling 

Before deriving the Inverse of A, Equation (5) Is modified so that the 
covariance matrix < z z"> has equal diagonal elements and the variables 
are scaled dimensionless quantities of order one [Wunsch, 1978]. The 
dimensionless weighting matrix W defined by 

a 2 At.j 2 

( 8 ) W i j = — y ; i , j = l,2n 

2n 

where 2 = J- z 0 ./& t. 2 , 6 u = 0 if i t j, and 6. . « 1 if i = j, is used 
to transform Equation (5) as follows: 


( 9) 





z 


where . = / W. . . Observe that now z z^ > = o 2 X- To avoid 

i j i j — — — 

1/2 

unnecessary notational complication, however, the factors jy will be 
henceforth dropped from Equation (11). 

1/2 

The purpose of the weighting matrix W is to render each component of each 

observed stake displacement rate in a form having the same uncertainty. In practice, 

6x . 6y . 

the displacement rate vectors ( » ) i = l,...,n are measured in polar 

6Rj 60. At.. 

coordinates ("Xt - ’ 7T - ) i=l,...,n where the uncertainty of 66. is considerably 
6R_- 1 1 

higher than — . This disparity results from the practical limitations of 

K i 

conventional theodolite and electronic distance measuring devices. Typically, 

_ 9 

accuracies of ± 0.5 x 10 m and ± 10" are attainable for 6R.. and 66^, respectively. 

Unless the outlying stakes are placed closer than -100m from the central stake, 
these disparate measurement uncertainties will render the uncertainty of the 
tangential component of ( fix^ , 6y 1 - ) greater than that of the radial component. 
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1/2 

The advantage of using the weighting matrix W is that it allows the data 

elements of d representing radial displacements to weigh more heavily than 

those representing tangential displacements. 

Nondimensionalization of variables is accomplished by defining the 

2^2 

length and strain-rate scales L, * = j] » and by transforming the 

variables as follows: 


(10) 


and , 


A. . 

TJ 

L) A. . 

1J 

9U i 

du. 

• i 

9Xj 

* C 9x j 

Ax- 

«y £ \ (L 

k 



(6 v 6 V 


2 


£ 


° Z £ 


9U. 


where A. . » 6x , <5y , and z are now of order one. Equation (5) may be 

expressed in terms of the non-dimensional uncertainty parameter / —2 | : 

\ e L At / 


(id 


A s 


L At 


Scaling Equation (9) in this way will simplify the forthcoming discussion 
on how to best select At and L from prior estimates of * and o. 

Singular value decomposition 

The estimate S that minimizes the error | | A - d || defined by: 
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Is obtained by the singular value decomposition method [Lawson and 
Hanson, 1974; Wunch, 1978]. According to this method, A is decomposed 
into the product of three associated matrices that are more readily 
invertable: 

A = g r R T 

(13) ~ 

The 4x4 matrix r is diagonal and is composed of the eigenvalues A.. , 1*1,4 
associated with the following linear equations 

(14) A T A r = x/ r t • 1.4 

and 

x £ 2 a * 1 = 1>4 

0 £ " 5 » . • . » 2n 

The nomalized eigenvectors of A A T and A T A, £ and r , respecti vely, comprise 
the columns of Q and the rows of respectively. The inverse matrices 
and Q~* are easily obtained by taking the transposes of R^ and Q. This 
simplicity results from the orthonormality of the eigenvectors and r . 

The inverse of A is* 

(16) = _1 = 5 L' 1 ? T 

*The complex algebraic manipulations required to perform singular value 
decomposition of A are easily done by computer programs now available on most 
university computer systems [Dongarra and others, 1979]. 
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where r”* Is the 4x4 diagonal matrix composed of diagonal elements 
* 1/Xi 1 * 1,4. Equation (9) is 
solved by applying A”* to the observed data d: 

r I?) Se = A" 1 d = R r" 1 S T d 

The expectation value of the error variance, <(se-s)(se-s) T >. . i=l,4, is 
given by: 


(18) 


<(Se-S)(Se-S) T > 1 . = < ^ r' 1 z 2 T ^ r' 1 R T >. 


[ k rR Wk] 


(c L At)' 

Equation (18) states that the error of each component of Se depends on the 
eigenvalues of £ and the eigenvectors r. , i=l,4. 

Optimal rosette design 


The task of designing the best rosette amounts to choosing the adjustable 
parameters At, L and (x^.y^) i=l,n that satisfy a design criterion given by 
Equation (18) restated as follows: 


(19) 


Max 

i 


k 


2 

o 


(c' L it ) 3 


< E 


> 


2 

where E is the desired accuracy level (expressed in units scaled by c ). 

Here, o 2 is a parameter combining information about the surveying equipment 
and natural strain-rate variation, and l is estimated or determined from 
prior data. In practice, it is best to choose L as long as possible (-2 km); 
therefore only n and At will be considered adjustable in the following discussion. 
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To derive the design criterion, the n-leg rosette will be 


Idealized as a system of n outlying stakes placed regularly on a circle at 
angular Intervals of ae * ~ . I" practice, however, allowances must be 
made for deviations from this Ideal configuration. Assuming that fcR and Ue 
are the Intrinsic radial and tangential uncertainties of observed stake 
displacements, the non-dimensional and pre-welghted version of the matrix A 


Is given by: 


(19) 




( 6R 2 cos 2 1ae+L 2 6e 2 s1n 2 iae ) ** 


(fiR^sin^Ue+L^fie^cos^lae) 15 


cos lae j e l and l£n 
sin iae j*2 and 1<n 
0 j*3,4 and 1<n 

0 j e l,2 and 1>n 
cos lae j*3 and 1>n 
sin lae j a 4 and 1>n 


where <j = 


2n i=i I 


(6R^cos^iae+L^6e ? sin^iAe) ,s 


(6R 2 $in 2 ia6+L 2 6e 2 cos 2 iae) 




The value of Max [e R 2 /X 2 ][o 2 /(eUt) 2 ] expressed as a percentage of e is 

1 k ik k 

presented as a function of n and at in Figure (3) for parameter values L=2 km 
and c=l(f 10 sec' 1 . Observe that this -error value varies as (at/fT)’ 1 . As shown 
in Figure (3), a choice of n=13 would yield 3% accuracy within a time span of 
4 weeks. 

Demonstration 

A stake network planted upstream of the Crary Ice Rise on the Ross Ice 
Shelf demonstrates the data-analysis procedure associated with the n-leg 
rosette. Figure (4) shows the disposition of four 3-leg rosettes (O.E.F.G) 
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planted in 1983 as part of a field project by Robert Bindschadler and 
others (in preparation) to measure strain rates upstream of an ice rise. 

The position of a much older 3-leg rosette planted in 1973 (Gg) falls 
within the area enclosed by the four new rosettes. This older rosette 
was re-surveyed in 1974 as part of the Ross Ice Shelf Geophysical and 
Glaciological Survey (RIGGS) [Thomas and others, in press]; and yielded 
an accurate strain-rate determination (±9%) with which data from the new 
rosettes may be compared. 

The four rosettes planted in 1983 were each surveyed and re-surveyed 
within a time interval spanning 4.2 days. This re-survey was conducted 
primarily to check survey equipment and to scan for possible zones of 
ultra-high deformation. Here, however, the re-survey data will be used to 
test the ability of a 12-leg rosette to yield an accurate measurement in 
a short time interval. The imaginary 12-leg rosette is synthesized from 
the four smaller rosettes by treating each of their central stakes as one. 
Actual strain-rate gradients expected within the field area will degrade 
the accuracy of the synthesized 12-leg rosettes because the smaller 3-leg 
rosettes are separated from each other by up to 15 km. Nevertheless, the 
synthesization provides an otherwise unavailable opportunity to demonstrate 
a 12-leg rosette measurement. 

Figure (5) shows the stake disposition for the synthesized 12-leg 
rosette, and indicates the observed stake velocities relative to the 
combined central stakes. 

for two stakes corresponding to rosette “0", 
tangential displacements were not measured. Angular displacements required 
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to observe the tangential components would have been below the accuracy 
level (±10") of the theodolites used in this project, so any unknown 
tangential displacements were arbitrarily assumed to be zero. 

The results of the 12-leg rosette analysis are compared with the known 
strain rates in Figure (5) and in Table (1). The two principle strain-rate 
components (denoted ejj and * 22 ) given by the 12-leg rosette differ from the known 
values by 9 % and 14%, respectively; and their orientations are displaced by 6°. The 
unsatisfactory results produced by the 3-leg rosette (E) closest to Gg, 
using the data analysis procedure outlined by Thomas and others [in press], 
are also shown in Figure (5) and Table (1) as a contrast to the results 
of the 12-leg rosette. 

Practical considerations 

It would be a mistake to adopt rosettes with large n without first 
considering the limitations imposed by logistic support, survey technology 
and the various benefits of a multi-year field program not emphasized in 
this study. It is, of course, always desirable to plant rosettes with 
large n. This desire conflicts, however, with the need to minimize the 
field-party workload, especially when aircraft-support limits the time a 
field party can spend at each field site. If field programs encompass a 
large number of field sites, the 3 or 4 leg rosette may be best because 
field operations would likely span several years regardless of whether 
individual sites were revisited in the same, or in a succeeding, field 
season. For field programs covering a small area, such as that conducted 
upstream of the Crary Ice Rise, rosettes with greater numbers of stakes 
may be preferable so that the entire field project can be accomplished in 


several weeks. 
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In view of the need to reduce the field-party workload, several design 
modifications to the n-leg rosette may prove useful. Figure (1) shows, 
for example, how planting k central stakes and m outlying stakes will 
achieve a high data redundancy while reducing the total number of stakes and 
the surveying workload. Surveying the m outlying stakes from k independent 
central locations yields k»m=n independent measurements of stake displacement. 
Only m+k<n stakes are required by this non-ideal n-leg rosette; and if 
appropriate measures are taken (e.g., assigning a separate corner-reflector 
prism to each outlying stake) the surveying can be accomplished within a 
reasonable time frame. 

Conclusion 

The stake-network designs that best measure surface strain rates and 
vorticity encompass a large degree of redundancy. The more stakes included 
in each network, the less individual stake displacement errors will effect 
the overall measurement. The advantage gained at the expense of the 
redundancy is the reduced time interval between survey and re-survey. 

Redundancy may also be a key element in the design of future field 
programs based on remote sensing techniques such as airborne photogrammetry. 
Brecher [1982], for example, has used photographs of natural features to 
identify the motion of the Byrd Glacier. These natural features are 
available in virtually unlimited number; so provide, in principle, a 
natural equivalent of a large-n rosette. A shortened time interval between 
initial and final aerial photography missions may reduce the camouflaging 
of natural features by drifting snow; and may, therefore, ultimately 
achieve greater accuracy by virtue of higher data redundancy. 
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FI6URE CAPTIONS 

Figure 1. The idealized n-leg rosette considered in this study for the 
purpose of error analysis consists of n outlying stakes surrounding a 
central stake (0) at equal angular intervals and with equal radial 
separation. Typical rosettes used on the Ross Ice Shelf had n=3 [Thomas 
and others, in press]. Also shown above is a modified n-leg rosette that 
uses m+k<n stakes. This rosette design achieves n-legs by having k = ^ central 
stakes, and may be more practical under actual field conditions. 

Figure 2. Uniaxial strain rate measured on the Ross Ice Shelf at station 

C- 16 (81°11 1 38”S , 189°30 ’O^’W) is plotted as a function of distance along 

a line bearing 312°T. The relative ice-shelf surface elevation along the 

line is also plotted. Open circles denote stake positions. The error 

-4 

level of the observed strain rate is approximately 1.0 x 10 /yr. 

Figure 3. The maximum error of the velocity gradients, expressed as 
percent of *, is plotted as a function of n and At. Scale factors used in 
the calculation of these graphs are * = 10" /sec and L = 2.0 x 10 -3 m. 

Notice that 3 % accuracy can be achieved by a 13-leg rosette in 28 days 
or by a 3-leg rosette in 84 days. Under some circumstances, it may be 
preferable to deploy rosettes having n>15 so that all measurements can be 
accomplished within one field season. This diagram, and other similar 
ones, can be used by field glaciologists to design their field program to 
maximize the data return for a given logistic effort. 
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Figure 4. Four 3-leg rosettes (0, E, F, G) planted upstream of Crary 
Ice Rise on the Ross Ice Shelf were re-surveyed after an average time 
period of 4.3 days. A fifth rosette (Gg) planted during RIGGS provides 
an accurate measurement of the actual strain rate to be used as a reference. 
In this study, the four 3-leg rosettes are combined into one 12-leg rosette 
to see how the subsequent strain rate measurement provided by the 4.2 day 
period compares with the actual strain rate. 

Figure 5. Stake positions of the 12-leg rosette are mapped along with the 
relative stake velocities. The principal axes of the strain rate measured 
using the 12-leg rosette are compared with those measured during RIGGS (Gg) 
[Thomas and others, in press]. Also shown are the principle axes resulting 
from the 3-leg rosette E using a 3.9 day interval between survey and re- 
survey. Uncertainty of the observed velocities is variable. The shaded 
ellipse represents, for example, the velocity uncertainty for the three most 
outlying stakes. The ellipse should be oriented perpendicular to the 
position vector. 



Table 1 . Principal strain-rates j 


ind e 


22I measured by various 


rosettes 


Year of 

Survey Rosette 
1973 G g 


At 

1 year 


111 

-10 -1 

1.72*0.15x10 sec 


22 

-1.57±0.03x10" ,0 sec _1 


Bearing 
of e 

1 1 

50° 



1984 

12-leg 

4.2 days 1.B7 

-1.79 

44° 1.77x10 

1984 

E 

3.9 days 2.67 

-2.74 

42° 


sec 
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